
function[mySound] = reverb(data, fs)

    reverbs = 16;
    delay = 0.025;
    
    
    mySound = data;
    for i=1:reverbs
        %shift this copy of the sound to the left
        sampleOffset = uint32( i * delay * fs);        
        thisReverb = zeros(length(data) + sampleOffset, 1);
        thisReverb(sampleOffset+1:end) = data * (1/i);
                
        %compose it with the output
        mySound = mySound + thisReverb(1:length(mySound));
    end
end
